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ABSTRACT 

This paper reports a comprehensive study on the gravitational wave (GW) background 
from compact binary coalescences. We consider in our calculations newly available 
observation-based neutron star and black hole mass distributions and complete an- 
alytical waveforms that include post-Newtonian amplitude corrections. Our results 
show that: (i) post-Newtonian effects cause a small reduction in the GW background 
signal; (ii) below 100 Hz the background depends primarily on the local coalescence 
rate ro and the average chirp mass and is independent of the chirp mass distribution; 
(iii) the effects of cosmic star formation rates and delay times between the forma- 
tion and merger of binaries are linear below 100 Hz and can be represented by a 
single parameter within a factor of ~ 2; (iv) a simple power law model of the en- 
ergy density parameter flGw(f) ~ f 2 ^ 3 up to 50-100 Hz is sufficient to be used as 
a search template for ground-based interferometers. In terms of detection prospects 
of this background signal, we show that: (i) detection (a signal-to- noise ratio of 3) 
within one year of observation by the Advanced LIGO detectors (Hl-Ll) requires a 
coalescence rate of ro = 3 (0.2) Mpc~ 3 Myr -1 for binary neutron stars (binary black 
holes); (ii) this limit on ro could be reduced 3-fold for two co-located and co-aligned 
detectors, whereas the currently proposed worldwide network of advanced instruments 
gives only ~ 30% improvement in dctcctability; (iii) the improved sensitivity of the 
planned Einstein Telescope allows not only confident detection of the background but 
also the high frequency components of the spectrum to be measured, possibly enabling 
rate evolutionary histories and mass distributions to be probed. Finally we show that 
sub-threshold binary neutron star merger events produce a strong foreground, which 
could be an issue for future terrestrial stochastic searches of primordial GWs. 
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1 INTRODUCTION 

Compact binary coalescences (CBC), of binary neutron stars 
(BNS), stellar mass binary black holes (BBH) and black 
hole- neutron stars (BH-NS), are the most promising source 
of gravitational waves (GWs) for ground-based interferome- 
ters such as LIGCQ and Virgqj. Although GW detections 
have not been recorded so far, a few tens of detections 
per year should beco me possible when a dvanced detectors 
come online in 2015 (|Abadie et al.ll201Ch . While individu- 
ally detectable CBC events are expected within distances 
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of hundreds of Mpc, the superposition of the gravitational 
radiation from these sources over cosmological v olumes can 
form a GW background (GWB; IPhinnevI lioOll ) . This sig- 
nal represents another interesti ng target for the up- coming 
advanced instruments (see, e.g., lRegimbaull201ll:l Zhu et al.l 
| 2011bl: iMarassi et al.|[2011bl : lRosadoll201ll : IWu et al ll2012l : 
Kowalska et al.ll2012l . for the most recent studies). 

A GWB is generally characterized by the dimension- 
less energy density parameter !]gw(/), which represents the 
present-day fractional energy density in GWs as a func- 
tion of frequency /. In general, assuming Newtonian en- 
ergy spectra and circular binary orbits for all sources, the 
CBC background can be described by a power law func- 
tion J1gw(/) = £l a f a , with a = 2/3 and an amplitude Q, a 
determined by system masses, coalescence rates and their 
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evolution over cosmic time. Such power law models have 
been widely used in sea rches for stochastic backgrounds 
using LIGO/Virgo data dAbbott et all 120091; lAbadie et all 
I2012T ), in mock data challenges dRegimbau et al.ll2O120 for 
the third-g eneration detector, the planned Einstein Tele- 
scope fET: IPunturo et al.ll201oT ) , and in parameter estima- 
tion of a stochastic background l|Mandic et al.|[2012f ). 

In this paper we investigate two issues of impor- 
tance for stochastic searches with ground-based interferom- 
eters. Firstly we refine the power law model for the CBC 
background by using complete analytical waveforms that 
include post-Newtonian (PN) amplitude corrections, and 
observation-based parameterized models of NS/BH mass 
distributions. The aim is to investigate what information 
can be extracted from and to provide a ready-to-use Hgw(/) 
model for CBC background searches. We secondly consider 
an additional motivation to study the properties of an as- 
trophysical GWB (AGWB) - the fact that it could act 
as a foreground masking the primordial GWBs from the 
very early Universe. As the spatial distribution of the in- 
dividual sources produces time series of varying GW am- 
plitudes, the strongest signals, which would be detected as 
single events, can be subtracted from the data. Therefore, 
as demonstrated for the BNS population u sing the proposed 
Big Bang Observer (|Cutler fe Harmsll2006ft . a detector with 
high enough sensitivity, could remove a foreground entirely 
by subtracting all the individually identified component sig- 
nals. We show in this work that there is a significant resid- 
ual foreground in the (1-500) Hz frequency range from sub- 
threshold BNS merger events. Such a foreground should be 
considered in future ground-based stochastic searches for 
primordial GWBs and other AGWBs. 

The organization of this paper is as follows. In section[2] 
we review the theoretical framework for calculating QGw(f) 
and other quantities of an AGWB used in the literature. 
We show the consistency among different approaches and 
confront the issue regarding the inclusion of an additional 
(1 + z) factor highlighted in a number of publications. We 
also present a practical power law model of AGWBs. In sec- 
tion0we extend this model to the general case for the three 
CBC populations by considering the effects of cosmic star 
formation rates (CSFRs) and delay times. Then using com- 
plete waveforms we calculate semi-analytically S1gw(/) of 
the CBC background. We describe in section [4] a Monte- 
Carlo approach to calculate !Jgw(/) which allow NS/BH 
mass distributions to be included and then show how the 
information of mass distributions is encoded in background 
energy spectra. In section [5] we use the improved models to 
carefully evaluate the detectability of the CBC background 
signal and further investigate the construction of £2gw(/) 
templates for future detectors. In section [6] we simulate the 
residual foreground noise for ET through the subtraction of 
the individually detectable events. In section [7] we discuss 
the unique time-frequency statistical properties of the CBC 
background and the possible implications with respect to 
detection. We present our conclusions in section [8] 



2 PROPERTIES OF AN AGWB 

In this section we summarize the broad range of formalisms 
presently used by different authors to calculate J1gw(/) 



of an AGWB. We start from Phinney's practical theorem 
|Phinnevll2"00ll ') and compare it with various versions given 
in the literature. 

Firstly recall that Qcw(f) is defined as the GW en- 
ergy density per logarithmic frequency interval at observed 
frequency /, divided by the critical energy density required 
to close the universe today p c = 3HqC 2 /StvG with Ho the 
Hubble constant. It is straig htforward to compute this di- 
mensionless function as fsee |Phinnevll200ll . for details): 
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f z =f(l+z) 



where N(z) is the spatial number density of GW events at 
redshift z; the factor (1 + z) accounts for redshifting of GW 
energy since emission; f T = f(l + z) is the GW frequency 
in the source frame and dEaw /dhxf T is source energy spec- 
trum. The limits of the integral over z are given by z m i n = 
max(0, fr n /f - 1) and z max = min(z„, / r max // - 1) with 
z* signifying the beginning of source formation and f" lln 
and / r max the minimum and maximal source rest-frame GW 
frequencies respectively. 

It is convenient to replace N(z) in equation |T} with the 
differential GW event rate dN /dz = N(z)c4-Kr z , where r z 
is the comoving distance related to the luminosity distance 
through dL — r z (l + z). We then obtain another version: 
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The quantity given by the above integration, with units 
of er gcm~ 2 Hz" 1 s" 1 , is c alled the spectral e nergy density 
(e.g.. iFerrari et al.lll999al: iMarassi et alj|2009h or the inte- 
grated flux (e.g., lRegimbaull201ll : IWu et al.l l2012h . Its di- 
mension shows that it can be related to the specific inten- 
sity by integrating the latter over the solid angle. The first 
two terms inside the integral give the locally measured en- 
ergy flux per unit freq uency (or simply fluence) emitted by 
a source at redshift z (|Flanagan fc HughesHl99"a ): 
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dSdf Airrl V df x 
while dN /dz can also be written as: 
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R(z) dV 
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with the comoving volume element dV/dz given by: 
dz H(z)' 



(3) 



(4) 



(5) 



where the Hubble parameter H(z) = Hq^J^Ia + fi m (l + z)' A 
and r z = J Q Z cdz' / 'H(z'), assuming a standard ACDM cos- 
mology with parameters H p = lOO/i-kms" 1 Mp c -1 , h = 0.7, 
O m = 0.27 and fi A = 0.73 (|jarosik et alj|201lh . 

In equation Q we define R(z) = roe(z) (see, e.g. 



ICoward et"aLll200l1 - |Howell et alj|2004h . which gives the rate 
density measured in cosmic time local to the event. The pa- 
rameter ro is the local rate density, usually used to estimate 
detection rates for different detectors, and e(z) is a dimen- 
sionless factor which models the source rate evolution over 
cosmic time. The later is usually associated with the CSFR 
for stellar catastrophic GW events. 
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The factor (1 + z) in equation Q converts R(z) to an 
earth time based qua ntity. The statement that suc h a factor 
does not exist given in lde Arauio fc Miranda! (|2005h does not 
change the calculation of Hgw(/) as the factor appears addi- 
tionally in their equation for dE/dSdf . This caveat has also 
appea r ed in othe r publ i cations , e .g.. iRegimbau fc Mandid 
(|2008h ; IZhu et al.l l|2010l , |2011al lbh: Howell et all (|201lh . We 
correct it with equations ((3| and Q since they provide phys- 
ically correct estimates of the corresponding quantities. 

Combining equations <j2j - (J5j yields the compact form: 
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(6) 

Alternatively J1gw(/) can be calculated through the 
single-source characteristic amplitude h c (f) = f{\h(f)\) 
with (\h(f)\) denoting the frequency-domain GW amplitude 
(in Hz -1 ) averaged over source orientations. In this case one 
can use the following relation to replace equation ©: 



dEGW _ 7TC 3 2/ rs 

dSdf ~ 2G cU; " 



(T) 



The average over all source orientations for an inspi- 
raling binary, a rotating NS o r a ringing BH is given by 
l|Satnvaprakash fc Schut3l2009h : 
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where l is the inclination angle of the characteristic direction 
of the source, determined by the orbital or spin angular mo- 
mentum, with respect to the line of sight. This shows that 
(\h(f)\) is smaller than that of an optimally oriented source 
(i.e., l = 0) by a factor of 4/5. 

The one-sided spectral density of a GWB, Sh(f), can 
be conveniently compared w ith detector sen sitivities and is 
related to fi G w(/) through (jMaggiordl200Cl ): 



^(/) = Iff /" 3 «gw(/). 



(9) 



Note that assuming an isotropic GWB, a factor of 1/5 should 
be included to account for the average detector response over 
all source locations in the sky, when the above equation is 
used directly to compare Sh(f) with noise power spectral 
densities of L-shaped interferometers. For instruments with 
non-perpendicular arms, it becomes sin 2 £/5 with £ being 
the opening angle between the two arms. 

Another important quantity of an AGWB is the (dimen- 
sionless) duty cycle, £, which describes the degree of overlap 
of individ ual signals in time domain. It can be computed as 
(see, e.g. JCoward fc Regimbaull2006t ): 
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where At is the average observed signal duration. A value of 
£ ^ 1 generally implies a continuous background. Here At 
is assumed to be frequency independent; If a dependence 
exists, the upper limit of the integration should be changed 
to z m ax. The above defined duty cycle may not be useful if 
there is significant frequency evolution of At, e.g., for CBC 
sources - this will be further discussed in section [7] 
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Figure 1. The dimensionless rate evolution factor e(z) based 
on different parameterize d models of CSFR. SF2 is taken from 
iPorciani fc Madaul J200l{^ and we use the " low r ate" gamma ray 
bursts derived model of [Robertson fc Elli 1 l|2012f l. 



2.1 A practical model 

We now derive a practical model for AGWBs formed by 
sources for which the gravitational energy spectrum can be 
approximated by a power law function of frequency. Such a 
case is of particular interest because Qgw(/) ~ f a is natu- 
rally obtained when dE G w/df T — Af"" 1 with A the overall 
amplitude. Then equation (O has a simple form within the 
frequency range / r min < / < / r max /(l + 
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pc H 

where we have defined a dimensionless function: 
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: dz. 



(11) 



(12) 



Note that we have changed the lower and upper limit of 
the integral from z m i n and z max to and z„ respectively, 
as we focus on the particular frequency range / r mm ^ / ^ 
1 + Z*) where the power law relation applies. 

We define e(z) = p*(z)/p„(0), where p,(z) is the CSFR 
density (in M© yr -1 Mpc~ 3 ). The assumption made here is 
that the GW event rate closely tracks the CSFR, e.g., in 
core collapse supernovae related mechanisms; Otherwise ef- 
fects of delay times should be included as we show in section 
13.11 for CBC events. Note that rn is equivalent to the pa- 
rameter A used in some studies to represent the fraction of 
stellar mass converted to GW source prog enitors (see, e.g., 
iRegimbau fc MandidkoOSl ; IWu et alj|2012l ). As estimates of 
ro do not normally rely on measurements of p*(0); rather 
they can be based on independent observations or theoreti- 
cal calculations, we choose to treat ro as a free parameter, 
independent on the CSFR models throughout the paper. 

In this work we consider five parameterized 



forms of p*(z) d erived 
Porciani & Madau 



from 



2001; Hopkins & Beacom 



observations (see 
20061: 



Fardal et al.l 120071 ; IWilkins et all 120081 : iRobertson fc Ellis! 
20121 , for details). The corresponding models of e(z) are 
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shown in Figure [T] We therefore set z* as the maximal 
redshift for which the CSFR model is app l icable : z* = 15 
for the recent study of iRobertson fc Ellia l|2012h which is 
derived from gamma ray burst observations and z* = 6 for 
the other four models. 

Figure [2] (upper panel) shows J(cv) calculated for a — 
[0, 5] using the five models of e(z). Since all the current pre- 
dictions of AGWBs in the frequency band of terrestrial de- 
tectors indicate that J7gw(/) increases from ab out 10 Hz 
to se veral hundreds Hz (see, e.g., Figure 6 of Re gimbaul 
l201ll ). the chosen range of a is adequate for most of pos- 
sible scenarios, e.g., a — 2/3 for inspiraling compact bina- 
ries as mentioned e arlier, a = 2 for NS r -mode instabilitie s 
jOwen et all 1 19981 ; iFerrari et all Il999bl: IZhu et al.l l2011aT). 
and a = 4 for magnetars jRegimbau fc de Freitas Pachecd 
l200fj ; iMarassi et all l2011al ). The five curves of J (a) are 
within a factor of 2 around the average, for which a least- 
square fit is log[J(a)] = 0.04a? 2 + 0.3a + 0.35. 

For a power law energy spectrum, the total GW energy 
emitted in the frequency range (/ r mln , / r max ) is AEcw = 
A J f a ~ 1 df . We further define a dimensionless function: 



K(a) = A(100Hz) c 



AEgv 



(13) 



to obtain a practical form: 
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The function K(a) obtained while arbitrarily setting 
/ r min =10 Hz and / r ma3t =1000 Hz is shown in the lower panel 
of Figure [2] The least-square fit of the average over the five 
models of e(z) is given by K(a) = (1.2 — 0.04a)/(a 2 — 1.1a + 
2.4). Note that the chosen values of (/ r min , / r max ) correspond 
to a frequency band where ground-based detectors have sig- 
nificant sensitivities; One can calculate K(a) for a specific 
type of source using equation (|13[) . Figure [2] implies that: a) 
as a increases, the high redshift (z > 4) sources contribute 
more to the background; b) effects of the CSFR introduce 
uncertainties in the overall amplitudes of Ugw(/) within a 
factor of about 2 for a < 3 and up to 5 for larger a. 

Combining equations {T5| and @ gives Su (/) in a con- 
venient form: 
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Similarly, a convenient relation between the duty cycle 
^ of an AGWB and At and To can be obtained by combin- 
ing equations @, and (| 10|) and averaging over the five 
models of e(z) shown in Figure [T] 



C = 0.2 



At 



To 



(16) 



lsecj \ 1 Mpc^Myr" 

Such a relation shows if there is a continuous GWB formed 
by one particular type of sources. 

Equations (|11|) - H15I ) represent our practical power law 
model for AGWBs. The power law relation holds for the fre- 
quency range [/ r mm , /! max /(l + 2*)], where effects due to rate 
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Figure 2. The dimensionless function J(ct) (upper panel) given 
by equation H21 and K(a) (lower panel) defined by equation d 1 3 1) 
calculated for the five models of e(z) shown in FigurefT] The bold 
line is a fit of average over the five models. 



evolutionary histories and cosmology are linear and the asso- 
ciated uncertainties should be reasonably within a factor of 
2. The model allows quick evaluation of the background sig- 
nal strength using approximate estimates of To and Ai?Gw- 
As our knowledge improves our model can be easily modi- 
fied to provide templates for future stochastic background 
searches. In the following sections we will develop a ready- 
to-use model for the CBC background by considering addi- 
tional issues that have not been included here. 



ANALYTICAL 



THE CBC BACKGROUND: 
APPROACHES 



In this section we extend the derivation in subsection [2711 to 
obtain models for CBC events analytically. 

3.1 A simple power law model 

Previous calculations of J1gw(/) for the CBC background 
have employed the Newto nian inspiral energy spectrum, 
with the exception of B BH (|Zhu et al.ll2011b1 ; lMarassi et al.l 
l2011bl ; IWu et al.l [20121 ). Following the previous derivation, 
we present here a simple power law model generalized for 
three CBC populations. 

In the Newtonian limit, the GW energy spectrum for 
an inspiralling circular b inary of comp onent masses mi and 
77i2 is given by (see, e.g.. lThornelll987l ): 



2/3,^5/3 



dfr 



-fr 



■1/3 



(17) 



where M c is the chirp mass defined as M c = Mrj 5 ^ 3 , with 
M = mi + 7712 the total mass and 77 — mim,2/M 2 the sym- 
metric mass ratio. Inserting this into equation (|6]) and com- 
bining the expression of p c gives (/ r min < / < / r max /(l + z»)): 
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where we have defined a dimensionless quantity: 



equation (|18[) becomes 



J 1/3 



Jo ^/nZT^jTTW 



-dz. 



(19) 



To determine the applicable frequency range of the above 
power law relation, one has / r mm well below 1 Hz and /™ ax 
given by the frequency at the last stable orbit (LSO) during 
inspiral /lso — 4400 Hz/M with M in units of M Q . 

The newly defined quantity J2/3 differs from J(2/3) as 
given in equation (|12[) in the definition of e(z): for CBC 
events, effects due to the delay time td between the forma- 
tion and the final merger of binaries should be taken into ac- 
count. By assuming compact binary formation closely tracks 
the cosmic star formation, we define e(z) — p» iC (.z)/p« !t: (0) 
by introducing a p*-related quantity: 



)dt d 



(20) 



where P(td) and f m i n denote the probability distribution for 
and minimum value of td respectively. The upper limit of 
the integral i» corresponds to z». For CBC events, P{td) fol- 
lows a l/t d forrrjf] as suggested by lates t population-synthesi s 
studies on compact binary evolution l|Dominik et al.|[2012l ). 
The parameters 2 and z/ are the redshifts when a GW event 
occurred and the system was initially formed respectively, 
with corresponding time coordinates t z and tf. In our fidu- 
cial cosmology, td is given by the lookback time between z 
and Zf, integrating dz' /[(l + z')H(z')] from z to Zf. The 
term dt / jdt z — (1 + z)/(l + Zf) is included to convert a rate 
at tf (p*(zf)) to one local to t z (p. jS (z)). 

It should be menti oned that our equatio n (|20H is equiv- 
alent to equation (2) of lReeimbau fe Hughes! (|2009|) by not- 
ing that the (1 + z) factor cancels with the one i n equation 
( H). Th e additional (1 + z) term in equation (9) of lZhu et alJ 
2011bh was an error, and lead to a factor of 2 underestimate 



of the BBH background signal. 

For the five considered CSFR models and for a mini- 
mum delay time t m i n in the range of 10 — 100 Myr, J2/3 is 
well constrained within (1.3 — 2.6). It is roughly a factor of 
2 smaller than J(2/3) as given by equation (|12[l and shown 
in Figure [2] where no time delay is assumed. Our selected 
range of t m j n i s large ly consistent with results presented in 
iDominik et al.l (|2012h : see their Figures 14-17 for details. We 
note, however, that in some extreme cases t m - m for BBHs 
could be much higher, e.g., 500 Myr. This does not change 
our results significant ly as we will show be l ow. F or the com- 
monly used CSFR of iHopkins fc Beacoml (|2006l ) . J2/3 as a 
function of t m i n (in Myr) can be expressed as: 



J 2 /3(tmin) = 3.67- 0.85 (t min )° 



(21) 



for lOMyr ^ t m i n ^ 500Myr; increasing t m i n from 100 Myr 
to 500 Myr reduces J2/3 from 1.85 to 1.3. 

Replacing the constants with their numerical values, 



3 We note that P(td) of type la supernovae was observationally 
found to be consistent with the l/td prediction s for progenitors of 
the GW induced binary white dwarf mergers lIGraur et alj|2011t 
iMaoz et al]|2012l) . 
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Here we have replaced M c in equation (|18|) with [M c } 
to account for a distribution of system masses - the consid- 
eration of (M c 5 ' 3 ) other than (M c ) 5 ' 3 is based on the fact 
that !Jgw(/) is an average over individual energy spectra 
characterized by M^ 3 . As the differences between the two 
quantities are very small (as we will show in Table [2j, we do 
not attempt to distinguish between them and will use the 
term average chirp mass. Note that the CBC background 
signal contains information about the physical chirp mass of 
sources, while single event detect ions normally measure th e 
redshifted chirp mass M c (l + z) (|Cutler fc Flanagan! 119941 ). 

We have reviewed calculations of Qgw(/) f° r the CBC 
background th rough a simple power law model. The model 
extends that of lPhinnevI (|200lf ) by considering different rate 
evolutionary histories and combining uncertainties associ- 
ated with CSFRs and delay times into a single parameter 
J2/3. In the next subsection we will introduce some addi- 
tional inputs to produce more accurate estimates. 



3.2 Beyond a simple power law 

We consider for the first time new information of two as- 
pects to refine previous estimates: 

1) Observation-based parameterized models of NS/BH mass 
distribution - through a Monte-Carlo simulation (section [4]), 
we will in subsection 14.11 investigate how the spectral shape 
of the background depends on the mass distributions; 

2) Up-to-date complete waveforms for populations of BNS, 
BBH and BH-NS systems - these will show how well a CBC 
background can be approximated by a simple power law 
model in the ground-based detector band. 

The main parameters are summarized in Tables [T] and for 
the interested readers we provide an overview below. Unless 
we otherwise specify, we will use the in formation contained 
in Ta ble [1] and [3 and the CSFR of IHopkins fc Beacoml 
(2006) in the following sections. 



3.2.1 Observational inputs 

We consider the parameterized models of NS /BH mass dis- 
tribution recently derived from observational mass measure- 
ments. For NSs that are observed in double NS systems (with 
one or two pulsars), high- precision mass m easurements are 
available (see Table 1 in I6zel et all |2012|. and references 
therein), indicating a very narrow distribution. Using the 
observ ational data for the 6 double NS systems, lOzel et aU 
(2012) found that the NS mass distribution can be well de- 
scribed by a Gaussian with a mean /1 = 1.33M0 and a stan- 
dard deviation a = O.O6M0. 

In contrast with the consensus on the narrowness of the 
NS mass distribution, the BH mass measurements are sub- 
ject to much larger uncertainties, leading to a greater range 
in inferred distribution. Utilizing the maximal amount of 
observational information available for 16 BHs in transient 
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low-mass X-ray binaries, IQzel et all l|201Cf ) concluded that 
the underlying mass distribution can be best described by a 
Gaussian with n — 7.8Mq and a — 1.2Mq. More recently, 
iFarr et al.l(|201ll ) considered a broad range of parameterized 
models, and using a Bayesian model selection analysis, they 
found a Gaussian and a power law distribution are preferred 
for low-mass X-ray binaries, whereas an exponential distri- 
bution and a two- Gaussian model are favored if 5 high-mass, 
wind -fed X-ray binary systems were included fsee lFarr et alj 
l201ll , for details). 

Unless stated explicitly in Table [1] our considered mass 
(in M Q ) interval for NS (BH) is [1, 2] ([4, 40]). Given the 
adopted models of distributions, it is highly unlikely to ob- 
tain masses outside these intervals. We note that the exis- 
tence of a "gap" between the maximum NS mass and the 
lower bound of observationall y inf erred BH masses h as been 
suggested in lOzel et al.l |2010l) and lFarr et al.l (|201lh . Such a 
"gap" can not be attributed to observational selection effects 
as concluded in the former paper. In this regard, terrestrial 
advanced GW detectors will be able to resolve this problem 
through precise measurements of NS/BH masses from tens 
up to hundreds of detections of CBC events. 

The BH spin distribution is highly uncertain - cur- 
rently there have been only about 10 stellar m ass BHs with 
(model dependent) spin est imates available (Mill er~et alj 
l2009l ; lMcClintock et al.ll201ll ) . Considering recent results on 
the d etermination of the extreme spin of the BH in Cygnus 
X-l jGou et al.ll201ll ; lFabian et alil2012T }. we assume a uni- 
form distribution with spin parameter x — S a /m 2 between 
-0.95 and 0.95, where S a is the spin angular momentum and 
m is the BH mass and positive or negative value of \ implies 
alignment or anti-alignment between component spin and 
orbital angular m omentum. As most NSs are observed to be 
weakly spinning ^Manchester et all I2005I - ) , and the fastest 
spinning NS in double puls ar systems, PSR J0 737-3039A, 
has a spin period of 22.70 ms l|Burgav et al .120031 ) and equiv- 
alently % ~ 0-05 (|Brown et alj|2012h . we neglect the spin of 
NSs in our analysis. 

Observational NS/BH mass measurements were also 
used as inputs or calibrati ons in the popula tion- 
synthesis simula t ions a dopted bv lMarassi et al] (|2011bh and 
iKowalska et all l|2012fi . Results of these studies are based 
on chirp mass distributions of some simulated populations 
of CBC sources. We assume in this paper that compo- 
nents of coalescing compact binaries follow the observational 
mass/spin distributions. Note that: a) for BNS, simulated 
chirp mass distribution presented in iDominik et al.l (|2012T ) 
is also very narrow and should give similar results to what 
we will obtain in the following sections; b) Our adopted BH 
mass/spin distributions only apply to BHs in X-ray binaries 
and may not be representative for BBH and BH-NS systems. 

3.2.2 Up-to-date analytical complete waveforms 

Following IZhu et all (|2011bh . we use the phenomenologi- 
cal inspiral-merger-ringdown waveforms fo r non- precessing 
spinning BBHs presented in lAiith et al.l (|201ll ). In this 
model the TaylorTl waveform is adopted for the inspiral 
phase, with 1.5 P N order amplitud e corrections to the New- 
tonian waveform l|Arun et al.ll2009T ). We note that the wave- 
form model is calibrated against numerical relativity simu- 
lations in the parameter range of mass ratios between 1 and 



Table 1. NS/BH mass (m) and spin (x) distribution. 



X 


m 


(m) 


Ref 


NS 


^(1.33,0.06) 


1.33 


(1) 




AT(7.8,1.2) 


7.8 


(2) 


BH [/(-0.95, 0.95) 


Power law 


7.35 




Exponential 


10 


(3) 




Two-Gaussian 


10 





Notes: All values of mass are in Mq. N(fi, a) implies a Gaussian 
distribution with a mean fj, and a standard deviation a\ U (a, b) is 
a uniform distribution between a and b. The upper/lower bound 
of BH spin corresponds to the recently determined extreme s pin of 
the BH in Cygnus X-l l|Gou et al.ll201ll ; lFabian et alj|2012h : pos- 
itive or negative \ implies alignment or anti-alignment between 
component spin and or bital angul a r mom ent um. References fo r 
mass distributions: (1) IQzel et al.l [120121 ); (2) I6zel et al.l [|2010r i 
— Gaussian BH mass, which is used as our fiducial model; (3) 
IFarr et alj J201lf) - for the other three models of BH mass distri- 
bution, and we use the median values of parameters given in the 



P(m) 



for 6 ^ m 23; Exponential 



paper: Power law - 

— P(m) ~ e m ' m °, with mo = 4.7 for m ^ 5.33; Two-Gaussian 

- N(7.5, 1.3) and iV(20.4, 4.4) with weights 0.8 and 0.2 respec- 
tively. The power law model in Ref (3) has a slightly lower mean 
fi = 7.35A/q due to the exclusion of one low-mass X-ray binary 
system in their analysis as compared to Ref (2). 

Table 2. Information about CBC populations used in this work. 





waveform 


ro 

(Mpc -3 Myr -1 ) 


(Myr) 


<M C 5/3 ) 
(< /3 ) 


BNS 


TaylorT4 


1 


20 


1.276 










4.948 


BH-NS 


IMR 


0.03 


30 


4.734 
5.795 
5.779 










24.22 


BBH 


IMR 


0.005 


50 


21.86 
35.40 
35.29 



Notes: IMR - the phenomenological inspiral-merger-ringdown 
waveform f or n on-precessing spinning BBHs presented in 
Aiith et al. I (l201lh ; we also use this model for BH-NS as an ap- 
proxi mation to the type-II spectr um found in numerical simula- 
tions jShibata fc Tamguchill201ll) . For BNS waveform wc adopt 
the TaylorT4 f o rmula with 3.0 PN amplitude accuracy given in 
iBlanchet et al. I feOOSl). Val ues o f ro correspond to the realistic 
estimates in lAbadie et al. I ll2010l) . imin given here is used as the 
fiducial valu e, based on the stand ard Submodel A for solar metal- 
licity Zq in Dominik et al. I l l2012l) - see Figure 8 therein; we also 
consider a range of 10-100 Myr to account for uncertainties. The 
quantities (M^ 3 ) are calculated using mass distributions pre- 
sented in Table [T] and assuming component masses are uncor- 
rected and follow the same distribution for BNS and BBH; four 
values for BH-NS and BBH are given in order from top to bottom 
as for a Gaussian, Power law, Exponential and Two-Gaussian BH 
mass distribution. We note that the quantity (M c ) 5 / 3 is smaller 
than (MrJ 3 ) by < 1% for BNS and the first two entries of BH-NS 
and BBH, and about 2% (4%) for the other BH-NS (BBH) values 
— we go with the latter quantity throughout the paper, but also 
use the former when comparing with other studies (in which case 
we neglect their differences). 
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4 and \ between -0.85 and 0.85, but we employ it for slightly 
broader parameter space. Our calculations can be improved 
once more accurate and general models become available. 

As no phenomenological complete waveforms are cur- 
rently available for BNqj and BH-NS systems, we consider 
analytical models that approximate the waveforms given 
by numerical relativity simulations. For BH-NS, we use the 
same model as that for BBH. The justification for our choice 
is two-fold. Firstly, the type-II spectrum found in numerical 
simul ations is similar to that of a BBH with the sam e mass 
ratio (|Shibata et al.ll2009l ; lshibata fc Taniguchil201l] ), show- 
ing a clear signature of inspiral, merger and ringdown. This 
happens primarily for larger mass ratios (> 3 — 5) when the 
smaller NS is simply swallowed by the BH. For the NS/BH 
mass distribution used in this work this condition is largely 
fulfilled. Secondly, PN amplitude corrections and effects of 
BH spins can be included by using the adopted BBH model. 

For BNS, we use the TaylorT4 point -particle wave- 
form with 3.0 PN order amplitude accuracy (Bl anchet et al.l 
2008). We apply the waveform up to 5000 Hz to ac- 
count for a realistic cutoff of the complete spectrum. Com- 
parisons between the TaylorT4 waveform and numerical 
relativity results generally indicate th at the former un- 
derestimates the po st-merger emission dKiuchi et al 1 l2009l : 
iFaber fc Rasioll2012l ). Recentlv lBauswein et all (|2012f l found 
that the generic outcome of two 1.35M0 NS mergers is the 
formation of a deformed differentially rotating massive NS, 
and that violent oscillations of the merger remnants lead to 
a pronounced peak in the GW spectra. We note that the 
peaks shown in this work are sharper than res ults obtained 
in full general relativistic simulations (see, e.g-. lKiuchi et all 
120091 ; iRezzolla et al.ll201()h . largely due to a different numer- 
ical treatment. 

Figure [3] shows the energy spectra for a BNS of equal 
mass 1.33M0. We consider a simple Gaussian spectrum to 
investigate the possible contribu tion from t he po st-merger 
emission to the GWB. Following Izhu et all (|2010l ). we take 
the form of dE G w/df = Aexp[-(f - / peak ) 2 /2A 2 ] where A 
arbitrarily set to be twice that of TaylorT4 waveform at f 000 
Hz, /peak = 1840 Hz and A = 250 H z - / De ak corresponds to 
the lowest value given in Table 2 of H auswem et al.l (|2012T ) 
and we use a much higher width A. The chosen parame- 
ters give a optimistic representation of post-merger emis- 
sion because: a) depending on NS equation of state, the 
peak frequency can be higher (up to about 4 kHz), together 
with narrower peaks, making it harder to detect (in terms of 
both single events and the contribution to a GWB); b) for 
the case of prompt BH formation (mainly for larger binary 
masses) the peaks are much smaller. Due to these uncertain- 
ties, the above mentioned Gaussian spectrum will be used 
only for semi-analytical calculations presented in the next 
subsection. 



3.3 Semi-analytical results 

Figure|4]compares three models of ilGw(/), using a NS (BH) 
mass of 1.33 (7.8) Mq and zero BH spin, and assuming that 
sources of each population have the same mass/spin values: 



4 Work is in preparation by the numerical relativity group at AEI 
and collaborators (L. Rezzolla 2012, private communication). 



TaylorT4 
Post-merger 

r 1/3 




10' 



10 



10 



/ (Hz) 



Figure 3. The GW energy spectrum for a BNS of equal mass 
1.33Mq calculated using the TaylorT4 formula. The post- merger 
signal is represented (optimistically) by a Gaussian spectrum cen- 
tered at around 2 kHz. Also shown arc curves of/" 1 / 6 and /- 1 / 3 
- / -1 / 6 corresponds to the Newtonian inspiral spectrum given by 
equation l|17[ l and / -1 / 3 shows the gradient at around 1 kHz. The 
curves are displayed as the square root of dEQyf/df in order to 
be directly comparable to the quantity h c g = f\h(f)\ commonly 
used in the numerical relativity community. 



1) A semi-analytical model calculated using equation ((6]l 
with complete waveforms described in subsection 13.2.21 

2) A Newtonian model based on equations © and (|17[) . and 
assuming / r max = /lso; 



3) A simple power law model based on equations (12 1 p and 
(J22J with an upper frequency cutoff /lso/5. Note that: a) 
an exact power law relation applies only for / /lso/(1 + 
z») with z, = 6; we empirically set the cutoff at /lso/5 
since the function inside of the integral in equation (|19|l has 
negligible values for z ^ 4; b) As mass distributions are not 
considered here, (M^ 3 ) in equation ([22} becomes M c 5/3 with 
M c determined by two component masses mentioned above. 

The following features can be observed from Figure [4] 
a) Newtonian models can be perfectly described by simple 
power law models up to /lso/5, about 300 Hz, 80 Hz and 
60 Hz for BNS, BH-NS and BBH respectively, as suggested 
in the previous paragraph; b) Semi-analytical models start 
to drop slightly below a / 2 / 3 power law from a few tens Hz 
due to PN amplitude corrections; c) Newtonian models give 
incorrect peaks and the followed abrupt decline because of 
the exclusion of post-inspiral emission. 

For BNS, we specifically show that: a) the power law 
index of H G w(/) drops from 2/3 (/ < 100 Hz) to 1/3 
before peaking at around 1-2 kHz; b) if the post-merger 
emission is included in the form of Gaussian spectra, the 
peak of Qgw(/) can be considerably enhanced while the 
low frequency part (< 300 Hz) stays at the same level. We 
will show, however, in section [5] that the contribution from 
post-merger emission to the background is unlikely to be 
detectable even with ET. 

Figure [5] shows S7gw(/) of the semi-analytical models 
using 10 different forms of e(z) based on the five CSFR mod- 
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- Semi-analytical 
-Newtonian 

.J2/3 




Figure 4. The energy density parameter S1gw(/) of the GWBs 
formed by three CBC populations (BNS, BH-NS, and BBH) cal- 
culated using complete waveforms (semi-analytical), compared 
with Newtonian models and simple power law models / 2//3 (see 
text). For BNS, the bump at around 1 kHz corresponds to the (op- 
timistic) contribution from the post-merger emission represented 
as a Gaussian spectrum shown in Figure [3] The BBH curves are 
scaled up by a factor of 4 to separate the three groups of curves, 
which is also the case in Figures 15171 



els and two minimum delay times t m i n = 10, 100 Myr. Two 
main results are: a) for each population, different curves fol- 
low the same gradient up to 100-200 Hz. This is in agreement 
with our derivation in subsection 13. II where we have shown 
within this frequency range effects of the CSFR and delay 
times are linear. There is a degeneracy between the CSFR 
and i m i n : to break this degeneracy a fully reconstructed e(z) 
and precise CSFR measurements are required; the former 
could become possible if we can efficiently detect most of 
the individual events out to high redshift, e.g., as we will 
show in section[5]for the BBH population; b) the only distin- 
g uishable feature comes fro m the adoption of CSFR model 
of iRobertson fc Ellis! (|2012l ); the relatively high CSFR from 
z — 3 up to 2 = 15 shifts the peaks of Ogw(/) to lower 
frequencies and suppresses the post-peak amplitudes. 



4 MONTE-CARLO SIMULATION 

In this section we describe a Monte-Carlo simulation ap- 
proach to calculate £7gw(/) of a general AGWB. This will 
allow us to investigate two important aspects of the CBC 
background in the next two sections: 

1) The dependency of f^Gw(/) on NS/BH mass distribu- 
tions; 

2) How much of the CBC background can be removed 
through single events detections to allow greater accessibil- 
ity to primordial GWs from the big bang. 

Combining equations ([2|-((5| and (0) yields: 



fiGw(/) = 



1 ro 2n c 




Figure 5. As in Figure [4] but only shows the semi-analytical 
models calculated for 10 different rate evolution models based on 
the five CSFRs shown in Figure [T] and two minimum delay times 
fmin = 10, 100 Myr. There is a degeneracy between CSFR model 
and t m i n below 100 Hz and the un ique signature at aroun d the 
peaks is due to the CSFR model in IRobertson fc EUid |2012h . 



where we have defined 



r 2 z e{z) 



(i + z)^n A + n m (i + z) :i ' 



(24) 



assuming that there is no correlation between the rate evo- 
lution and source intrinsic parameters. Here we use an "ob- 
served" frequency domain amplitude in which detec- 
tor responses are included, which will enable us to estimate 
the signal-to-noise ratios of individual events in section [6] 

The discrete version of the integration in equation 
is a sum over events distributed in redshift, leading to: 



fiow(/) 



1 r 2tt 2 c 3 



H 

JV, 



G 



(25) 



1 JV >>»c 



! ffQi) 



f\h(f)\ 2 g(z)dz, (23) 



where i denotes the i-th event; Qi contains the intrinsic 
source parameters, which in our case includes binary com- 
ponent masses m\ and m\, and BH spin parameters x\ an d 
X\\ the parameter /, is the observed frequency, which de- 
pends on both Qi and Z;; the function P(z) is the proba- 
bility distribution function of source redshift z; and 7V mc is 
the number of events in our Monte-Carlo simulation, cho- 
sen to be 10 6 - the approximate expected number of BNS 
merger events within z* in one-year observation. Note that: 

a) f2Gw(/) does not depend on JV mc or an observation time; 

b) l/[Nm C P(z)] plays the role of dz in the integration of 
equation (123 p - 1/P(z) is essentially a weight used when 
calculating the average over individual sources; it simply 
becomes the length of integration (z max — z m in) without any 
prior knowledge of source redshift distribution, e.g., for semi- 
analytical integrations of section|3]and other similar studies. 

The effective distance D c ff of the source is related to its 
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luminosity distance through (|Allen et alj|2012h : 



D eg = dh 



., / 1+ COS 2 A 2 2 2 

F + I ] + F x cos t 



-1/2 



(26) 



where _F + and F x are the antenna pattern functions for + 
and x polarized GWs respectively, depending on source po- 
sition with respect to the detector (described by the right 
ascension 9 and declination <f> of the source) and the po- 
larization angle ip; t is the inclination angle. We use the 
forms of F+ and -Fx for a L-shaped interferometer like 
LIGO; when averaged over uni formly distrib uted 8, <j> and if), 
{Fl) = (Fx) = 1/5 (see, e.g.. lThorndll987l ). Therefore, re- 
sults obtained through equation (|25|) which will be presented 
in the following sections are a factor of 5 smaller than those 
calculated in section @. 

An initial step in performing a Monte-Carlo simulation 
is to construct probability distribution functions of param- 
eters. Here we use the NS/BH mass and BH spin distribu- 
tions given in Table [T] and further assume that mi and rti2, 
Xi and X2 are uncorrelated. The parameters cos 6, 4>/n, ipf-K 
and cos i are all uncorrelated and uniformly distributed over 
[—1,1], where the consideration with cos# and cost ensures 
that individual sources and the direction of their orbital an- 
gular momentum are uniformly distributed on a spherical 
surface. The function P(z) is obtained by normalizing the 
differential event rate given in equation Q. The parameter 
d,L is dependent on redshift z, given by r z (l + z). 

Our final results are obtained using an average of 10 in- 
dependent realisations of the Monte-Carlo simulation. Nu- 
merical error in our simulation, defined as the relative varia- 
tion (between each realisation and the average) of reference 
values of J7gw(/) at 100 Hz, are within 3%. As individual 
sources contribute to the background through a / 2//3 power 
law below 100 Hz, the outputs of different realisations vary 
by a small linear scaling in magnitudes of ilaw(f)- Note 
that our results represent the average background energy 
spectra. The actual background signal can vary considerably 
from the average depending on the time-frequency proper- 
ties as we will discuss in section [7) 

We shall present our results as both: a) a full back- 
ground - calculated using equation Q25p for each population, 
which is the background signal we will be searching for; b) 
a residual foreground - in the summation of equation (|25|) 
individual events above the detection threshold for ET are 
discarded, which represents the residual noise due to sub- 
threshold sources. Result a) will be presented in section |4~T1 
and a) and b) are compared in section [5] with a focus on b). 
For completeness, we also present an example of a simulated 
time series due to the BNS population in section [7] 



4.1 Backgrounds encoded with mass distributions 

We compare the numerical results of S1gw(/) for each CBC 
population, with the semi-analytical model presented in sub- 
section [3]3] and a simple power law model described in sub- 
section [371] The function Q,aw(f) for the later two models 
are scaled down by a factor of 5 due to an average of detec- 
tor response over source locations, to be comparable with 
the numerical results. 

Figure [6] shows such a comparison in the case of Gaus- 
sian mass distributions. The very narrow distribution of NS 




Figure 6. The CBC background energy spectra calculated nu- 
merically assuming Gaussian mass distributions are compared 
with the semi-analytical models and simple power law models 
/ 2 / 3 (up to 200 Hz) shown in Figure [4] Curves of the later two 
models are scaled down by a factor of 5 due to an average of detec- 
tor response over source locations, which means !Jqw(/) shown 
in this plot (and also in Figures [7] and lilt are 5 times smaller 
than those calculated in section [3] 



- Gaussian 

- Power law 
Exponential 

- Two-Gaussian 




Figure 7. As in Figure[6] but showing only the numerical results 
of BH-NS and BBH for 4 models of BH mass distribution. Curves 
are scaled linearly to the same (M^ 3 ) as that of a Gaussian 
distribution (see Table [S). 



masses has negligible influence on the BNS background - the 
numerical result perfectly matches the semi-analytical model 
except a slightly broader shape around the peak, while ef- 
fects of a Gaussian BH mass distribution are moderately no- 
ticeable for BH-NS and BBH - the reduction from the / 2/3 
curve is partly alleviated due to the contribution from the 
merger-ringdown emission of more massive systems. Unless 
otherwise stated we adopt the numerical models assuming 
Gaussian mass distributions in the following sections. 
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Figure [7] illustrates the numerical models of BH-NS and 
BBH for different BH mass distributions. Curves are scaled 
according to the individual values of (Aff 3 ) as given in Ta- 
ble [2] for each distribution to ensure that all have the same 
value as that of a Gaussian distribution. Two groups are 
clearly distinguishable above 200 Hz - one containing low 
mass BHs only (Gaussian and power law) and another for a 
broader distribution when high mass BHs are also included 
- different peak frequencies are due to variations of average 
total masses and slightly distinct spectral width comes from 
the different degree of concentration of the distribution. The 
amplitudes of S1gw(/) up to 100 Hz are very similar to each 
other (within the 1% numerical errors), agreeing with the 
expected dependency on (Mc^ z ). 

Two main conclusions from Figures [6] and [7| are: a) mass 
distribution plays a role only through {M^ A ) in the low- 
frequency (< 100 Hz) power law part o f the backgroun d 



energy spectrum. This was mention ed in Wu et al. J2012h . 



and also confirmed independently in iKowalska et all (12012T ) 
where a power law relation was obtained using mass distri- 
butions derived from population-synthesis simulations; b) it 
could become possible to probe mass distributions through 
stochastic background measurements, e.g., peaks shown in 
Figure [7] (once measured) will provide information about 
the average total mass and the degree of concentration of 
the distribution. 

In our calculations we do not consider more sophisti- 
cated distributions of BH spin and mass ratio. These two 
parameters play a minor role (compared with M c ) above 
a few tens Hz in our adopted waveforms. Therefore, our 
results will not be affected significantly as long as their 
true distributions are not highly asymmetrical. An addi- 
tional effect due to orbital eccentricity is not relevant as 
the orbits of coalescing compact objects are expected to 
circularize before t heir GW signals enter the gro und-based 
frequency window l|Brown fc Zimmerman! I2OIO1 ). We note 
that dynamically formed BBHs, of which the population 
is not considered in the current work, may be highly ec- 
centric and could merge before their orbits are c ircularized 
1 Benacquista fc Downingll201ll ). As mentioned in lZhu et al.l 
I 2011bl ). such a population, possibly with much higher av- 
erage masses, could provide considerable contribution to a 
GWB. 



5 ISSUES ON THE DETECTABILITY 

In this section we first update previous estimates on the 
detectability of the CBC background for second and third 
generation terrestrial detectors, using improved background 
models. By considering practical issues in detection and pa- 
rameter estimation, we further discuss the choice of fiGw(/) 

templates for d ata analysis . 

Following IZhu et al.l (|2011br ), we consider five 
advanced detectors - advanced LIGO (aLIGO; 
lHarrv fc LIGO Scientific Collaboration! |201Ch at Han- 
ford (H) and Livingsto n (L), advanced Virgqj (V) in Italy 
jAcernese et al.l 120091). KAGRA (K, previous l y kn own 
as LCGT. iKuroda fc the LCGT Collaboration! I2010I ) in 



Japan, and AIGO in Australia (Ba rriga et al.l l201Cf ). as 
well as ET0 for which we conside r two configurat ions - 
ET-B (|Hild et al.ll200g| ) and ET-D (|Hild et al.ll201ll ). Note 
that the inclusion of IndlGCQ in India should have similar 
contribution as AIGO. Unless otherwise stated, we use the 
sensitivity curves of aLIGO for the zero detuning, high 
laser power configuration (see the public LIGO document 
T0900288 for detail©, and of KAGRA for the broadband 
configuration Q. We assume that AIGO has the same 
sensitivity as aLIGO. The target sensitivities of these 
detectors are shown later in Figure [12] 

5.1 Signal-to- noise ratios 

The optimum detection strategy for a stochastic GWB is 
to c ross-correlate the outputs of two or mor e detectors (see, 
e.g.. lAllen fc Romanoll 19991 ; lMaggiorell2000l '). Strictly speak- 
ing, the CBC background is not a stochastic background 
in the sense that individual signals do not s ufficiently over - 
lap in time-frequency space, as suggested in iRosadol (|201 lh 
and we will also discuss in section [7] Nevertheless it ha s 
been shown, both theoretically (|Drasco fc Flanagan! I2003T ) 
and experimen tally (through data a nalysis exprement on 
simulated data; [Regimbau et al.1120121 ) . that the cross corre- 
lation method works nearly optimally in the non-Gaussian 
regime, because through long time integration it is always 
possible to obtain a sufficiently large number of signals in 
a frequency interval and "form" a Gaussian background for 
which the cross correlation statistic applies. Therefore, we 
consider this standard method to assess the detectability of 
the CBC background for future detectors. 

The optimal signal-to-noise ratio (S/N) obtainable by 
t wo-detector cross correla tion is given by (e.g., equation 3.75 
in lAllen fc Romano! fl999l ) : 



(S/N) 2 



2T 



df 



(27) 



http: //wwwcascina. virgo.infn.it/advirgo/ 



S nl (f)Sn2(f) 

where 7(/) is the normalized overlap reduction function, 
which accounts for the sensitivity loss due to t he separa- 
tion and relative orientation of the two detectors (jFlanaganl 
1993). For co-located and co-aligned detectors, 7(/) = 1. 
The one-sided noise power spectral densities of the two de- 
tectors are given by S n i(f) and S„2{f), and T is the in- 
tegration time (set to be one year). Note that: a) we have 
substituted J7gw(/) with the spectral density Sh{,f) through 
equation ((9)1 to obtain a more intuitive format; b) the factor 
1/5 for Sh(f) is ignored, as our numerical results of Qgw(/) 
includes effects of detector response. We us e the y(f ) for the 
10 pai rs of advanced detectors presented in lNishizawa et all 
(2009) and adopt the form of ET fo r two V-shaped detect ors 
separated by 120° (see Figure 8 of lRegimbau et ai1l2012l ). 

As we will be observing a GWB due to all possible con- 
tributions of CBC sources, we calculate the S/N of the total 
background from the three CBC populations without con- 
sidering other types of sources. The background spectrum 
of the total CBC background is simply the sum of that of 
each population, as shown in Figure PTT1 



6 http:/ /www.et-gw.eu/etsensitivities 

7 http://www.gw-indigo.org 

8 https:/ /dcc.ligo.org/cgi-bin/DocDB/ShowDocument?docid=2974 

9 http:/ /gwcenter.icrr.u-tokyo.ac.jp/en/researcher/parameter 
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Table 3. S/N of the CBC background for advanced detectors. 



Pair 


A-H 


A-K 


A-L 


A-V 


H-K 


a 


0.43 


0.11 


0.46 


0.10 


0.06 


b 


(0.43) 


(0.17) 


(0.46) 


(0.18) 


(0.10) 


Pair 


H-L 


H-V 


K-L 


K-V 


L-V 


a 


1.05 


0.12 


0.02 


0.12 


0.14 


b 


(1.05) 


(0.23) 


(0.03) 


(0.29) 


(0.25) 


Case 


HH 


KK 


VV 


Ca 


Cb 




2.76 


1.82 


1.31 


1.26 


1.34 



Notes: Assuming that the background is contributed by the three 
populations of CBC sources (BNS, BH-NS, and BBH), we present 
S/N of cross correlating different interferometer pairs from a 
worldwide network. Results are shown for: a - individual detector 
sensitivities (but note that A, H and L have the same aLIGO sen- 
sitivity); b - assuming all detectors have aLIGO sensitivity. The 
motivation for calculating b is to investigate effects of the overlap 
reduction function 7(/); HH, KK, and VV assume 7(/) = 1 for 
aLIGO, KAGRA and advanced Virgo respectively; Values of Ca 
and Cb are just the square roots of the quadratic sum of the 10 
values in a and b respectively, which can in principle be achieved 
by optimally combining measurements from multiple pairs of de- 
tectors. The improvement of Ca and Cb on H-L is not appreciable 
due to the suppressing effects of 7(/). 

We present in Table [3] values of S/N calculated for ad- 
vanced detectors. We consider three cases: 

1) Cross correlation between pairs of advanced detectors us- 
ing real 7(/) and the individual sensitivities of each detector; 
we additionally assume all detectors have the same sensitiv- 
ity as aLIGO to evaluate the effect of 7(/); 

2) Assuming 7(/) = 1 for pairs of detectors with the sensi- 
tivity of either aLIGO, or KAGRA or advanced Virgo; 

3) An optimal combination of cross correlation statistics for 
10 pairs of advanced detectors, for which (S/N) 2 is simply 
the su m of those calculated i n case 1) (see, e.g., equation 
5.46 in lAllen fc Romano|[l999h . 

Note that case 3) is mathematically simple but requires 5 
advanced detectors to be simultaneously online. 

Our results show that: a) among the 10 pairs, H-L per- 
forms the best in terms of detecting a CBC background, 
giving a S/N of 1, while the lowest value of S/N is only 0.02. 
Assuming the aLIGO sensitivity for all detectors only in- 
creases the lowest value to 0.03; b) the improvement from 
combining the network of advanced detectors is only ~ 30% 
on the best performing pair H-L, while assuming 7(/) = 1 
for aLIGO increases the S/N by nearly 3 folds. This is well 
below the ex pectation that t hese two should give similar 
improvement (|Wu et alj|2012l ; iKowalska et alj|2012r ). Such a 
pessimistic prospect is mainly due to effe cts of "/(f). This h as 
been pointed out in our previous studies |Zhu et al 1l2011al fbTl 
and will be discussed in more details below. 

The property of 7(/) is mainly described by its char- 
acteristic frequency / c har, given by / c h ar = C /(2|AX|) with 
|AX| the distance between two detectors, above which 7(/) 
decays rapidly towards zero. Among the 10 pairs of detec- 
tors, H-L has the smallest separation (|AX| = 3000 km), 
resulting in the highest / c har of 50 Hz, whi le values of / c har 
for o ther pairs vary from 10 Hz to 20 Hz (|Nishizawa et all 
2009). This, combined with the fact that advanced detectors 
have a low frequency seismic wall at about 10 Hz, can eas- 



ily explain the very small 30% improvement. Note that the 
overlap reduction function also depends on the relative ori- 
en tation of the two detect ors, and we refer interested readers 
to iNishizawa et all (2009) for discussion about the optimal 
configurations of (geographically separated) detector pairs. 

For ET, the CBC background can be easily detected, 
with S/N of 178 (350), 19 (38), and 15 (30) assuming ET-B 
(ET-D) sensitivity for the BNS, BH-NS and BBH population 
respectively - the factor of 2 increment from ET-D is due to 
greater sensitivity at frequencies below ~ 20 Hz. This im- 
plies that the detection prospects benefit significantly from 
improvement of low-frequency sensitivities. 

We note that the (optimistic) BNS post-merger contri- 
bution to the GWB, as shown in Figure [4] results in a S/N 
of only 0.43 (0.46) for ET-B (ET-D), implying that detect- 
ing the imprint of BNS post-merger emission on a GWB 
requires a coalescence rate at least 5 times higher than the 
realistic value adopted in this study. On the other hand, we 
quantify PN effects with the difference in S/N between a 
simple power law model and the numerical model shown in 
Figure IH] (both are cutoff at 200 Hz). We find that PN am- 
plitude corrections cause a reduction of S /N in the range of 
(5.5% - 8.5%) and (2.1% - 3.2%) for ET-B and ET-D re- 
spectively. ET-D is less sensitive to these effects as its best 
sensitivity is more concentrated at lower frequencies. 

5.2 Detection prospects for advanced detectors 

As the operation of advanced detectors is only a few year 
away, it is now important to carefully assess the detection 
prospects of the CBC background, which represents one of 
the most (if not the most) promising background sources. 
We look at this issue in much more depth by considering 
variations in both source parameters and detector configu- 
rations. 

For each population, Qgw{/) scales linearly with ro, 
for which the uncertainties are generally of orders of mag- 
nitude - much larger than those of other parameters. To 
determine what possible combinations of ro the total CBC 
background will be accessible to advanced instruments, we 
simply scale the numerical models for different values of ro, 
keeping all other parameters fixed. The considered range of 
ro (in Mpc -3 Myr -1 ) is 0.1 — 10 (pessimistic to optimistic) 
for BNS, and 0.005 - 0.30 (realistic to optimistic) for BBH, 
while the BH-NS rate is set to be the real istic 0.03 (all values 
taken from Table 4 of lAbadie et al.ll2010l ). 

The motivation of our choice regarding BBH and BH- 
NS is two- fold: a) it w as recently predicted, through popula- 
tion s ynthesis studies l|Belczvnski et af . 2010; Dominik et al.l 
2012) and empirical esti mation based on t wo observed BH- 
Wolf-Rayet star systems (|Bulik et al.ll201lT ). that r for BBH 
can plausibly be at the optimistic value adopted above; 
b) while the same population synthesis studies gave sim- 
ilar realistic rates o f BH-NS (see, e.g., Tables 2 and 3 in 
iDominik et~ai1l2012t ). a negligible coalescence rate for BH- 
NS was recently empirically determined by fo l lowin g the fu- 
ture evolution of Cyg X-l (Belczvnski et all 120111 ). In the 
current analysis for advanced detectors, the contribution of 
the BH-NS population is nearly negligible at the chosen rate. 

Figure[S]shows the detectable "rate space" for advanced 
detectors: a S/N threshold of 3 is used to indicate detec- 
tions, which corresponds to 95% detection rate and 5% 
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aLIGO, H-L 
advanced network 
aLIGO, 7 (/) = 1 



(|27|) we have a simple relation: 

3 

S/N = C 2/3 ^> fc <Af, 



5/3 



\kjk 

> J-. 



2/3 



(28) 




BNS r 

Figure 8. The detectable "rate space" of a CBC background, 
assumed to be contributed by three populations (BNS, BH-NS 
and BBH), for ground-based advanced detectors. The local coa- 
lescence rate ro is in Mpc- 3 Myr _1 with r for BH-NS fixed at 
0.03 (see text). A signal-to-noise ratio greater than 3 for one year 
observation can be obtained for rates above each curve assuming: 
a) cross correlation with aLIGO H-L; b) an ideal case of r f(f) = 1 
for two aLIGO detectors; c) an optimal combination of measure- 
ments of 10 pairs advanced detectors (advanced network). Note 
that this background will be easily detectable for ET, i.e., the 
corresponding curve is well below the origin of the two axes. 



false alarm rate. Note that: a) we have taken the integer 
3 for convenience, while the accurate number is 3.29, given 
by V2[eiic' 1 {2/3) - erfc~ 1 (2<;)] assuming a false alarm rate 
/3 = 5% and a d etection rate <; = 95%. Here /3 and ? need 
not sum to 1 (see I Allen &: Rom ano 1999], for details); b) one 
can lower the threshold on S/N with the cost of increasing 
/3, or decreasing <; or both. For example, a threshold of 2.56 
may be used if one chooses /3 — 10% and <; = 90%. 

In the c ase of 7(.f ) = 1 fo r aLIGO, our results are con- 
sistent with IWu et all (j2012T ). implying that the BNS pop- 
ulation alone may produce a detectable background signal 
at the realistic coalescence rate for two co-located and co- 
aligned aLIGO detectors. However, an important point here 
is that for a worldwide network of advanced detectors, the 
requirement for a detection is more than twice stronger. This 
motivates us to consider all three populations as a whole as 
they will be observed in reality. Figure [8] shows that some 
combinations of the BNS population and the BBH popula- 
tion can form a detectable GWB, while both of their indi- 
vidual contributions alone are not sufficient for detection. 
In practice, if a Qgw(/) ~ f 2 ' 3 power law background has 
been detected, one would certainly wish to determine the 
relative contribution from every possible population. 

The possible variation in (M c 5/3 ) (for BBH and BH-NS 
systems) and the effects of CSFR and delay times (which 
can be represented using the parameter J2/3) are not con- 
sidered in Figure [8] Combining the simple power law model 
given by equation (1221) , which will be shown to be a good ap- 
proximation in the next subsection, with equations ([9]) and 



where the indices k = 1,2,3 denote the three CBC popu- 
lations and C2/3 is a constant depending only on detector 
sensitivity and j(f) for different detector pairs used in cross 
correlation. For convenience we have omitted the division 
by the corresponding reference values as in equation ()22[) 
for the three parameters. Note that one needs to set the 
upper frequency limit of the integration in equation (|27[) at 
100 Hz so that the above equation is representative of re- 
sults obtained using numerical models of fiGw(/)- In fact, 
our Figure [8] can be easily reproduced by using equation 
1)281) together with values of {Ale' 3 } and t m i n (to obtain the 
parameter J2/3 through equation 1)211) ') given in Tabled 

Varying the values of (M c 5/3 ) for BBH and BH-NS 
from those of a Gaussian distribution (as assumed in Fig- 
ure [H]) to the highest entries in Table [2] increases the to- 
tal S/N by 5% (40%) assuming a BBH coalescence rate of 
ro = 0.005 (0.3) Mpc~ 3 Myr~ 1 and realistic values of ro for 
both BNS and BH-NS. Such a increment of S/N is smaller 
than 40% for higher coalescence rates of BNS and BH-NS. 
Therefore, our Figure[8]does not change appreciably for vari- 
ations of {Ml' 3 ) given the current observational BH mass 
estimates. Meanwhile, effects of CSFR and delay times could 
moderately degrade (i.e., no more than a factor of 2) the de- 
tection prospects, as our current choice gives J2/3 = 2.3 for 
the dominant BNS population, which is close to the high 
end of the range (1.3 — 2.6) obtained in section I37T1 

Note that for a putative population of dynami- 
cally formed B BHs in dense stellar clusters (see, e.g., 
ISadowski et al.ll2008l) or for the same field population (as 
considered in this study) b ut assumed to be fo rmed in 
low metallicity environments (|Dominik et aL|[2012r ), a much 
larger average chirp mass (M c ) up to about ~ 20Mq (in 
comparison to (M c ) ~ 7Mq for the mass distribution used 
in this work) was suggested to be possible. To allow for 
these possibilities, a plot of detectable ro — (M c ) space is 
useful. Such illustrative studies, which only apply to a sin- 
gle po pulation, have been presented for the BBH popula- 
tion in IZhu et all feoilbl) and for each of the three CBC 
populations in Iwu et all l|2012l ') - for a given S/N thresh- 
old, the scaling relation ro ~ {M c )~~ 5 ' 3 was shown to be a 
good approximation^"!. At a coalescence rate of the order 
10 _3 Mpc- 3 Myr- 1 (see, e.g.. iMiller fc Lauburd I2009T ). dy- 
namical formation scenarios should have similar contribu- 
tions to a GWB as the field population of BBHs considered 
in this study (at the realistic rate) and thus will not improve 
considerably the detection prospects for advanced detectors. 
Such a back-of-the-envelope argument also applies to bina- 
ries involved with one or two intermediate mass BHs - much 
lower rates cancel o ut the advantage of higher masses (see, 
e.g., Tables 8-10 of lAbadie et alll201Ch . Despite of the in- 
volved uncertainties, these systems should be a more inter- 
esting source for single event searches/detections (at least for 



10 We no t e that curves for advanced detectors in Figure 5 of 
IZhu et al. I feOllbh underestimate the detectability by a factor of 
4 due to the use of an old version of aLIGO sensitivity and one 
additional factor of (1 + z) in the calculation of !)qw(/). 
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advanced detectors), from which their very existence will be 
tested or the associated coalescence rates can be stringently 
constrained. 

Looking forward to the advanced detector era, it is 
now crucial to investigate how the detection prospects 
of the CBC background (which could be the first 
to be detected) can be enhanced. For co-located de- 
tectors, techniques to remove correlated environmental 
and/or instrumental noises will be required, and are 
current ly being developed in the LIGO /Virgo collabo - 
ration (|Fotopoulos fe LIGO Scientific Collaboration] [2008). 
Provided that no co-located instruments are available, de- 
tection of a GWB from CBC events will require higher co- 
alescence rates than what are presently thought to be re- 
alistic, i.e., r = 3(0.2)Mpc- 3 Myr~ 1 for BNS (BBH), or 
alternatively given the realistic rates an integration time of 
4 yea rs to obtain a S/N of 2, which was assumed as a thresh- 
old in lWu et al.l l|2012h . We note that a single-detector auto- 
correlation approach was recently proposed to be compara- 
ble in S/N to what is achievable by cro ss correlation of two 
co-lo cated and co-aligned detectors (|Tinto fc Armstrond 
2012), of which the feasibility needs to be further tested 
in realistic data analysis experiments. 

The above results have assumed standard versions of 
design sensitivities for advanced detectors. In practice, de- 
tectors can be tuned to different configurations for various 
purposes, e.g., to allow optimization for different searches. 
As the aLIGO H-L pair gives the majority of contribution 
to the network S/N for a CBC background, we consider here 
four additional tuning options of aLIGO (data for the corre- 
sponding sensitivities are available publicly at the link given 
in the beginning of this section, and we refer interested read- 
ers to the LIGO document T0900288 therein for descriptions 
and technical details): a) Zero-detuning, low laser power; 
b) Optimal NSNS, which is optimized to the BNS inspi- 
ral search; c) Optimal BHBH, which is optimized for 30-30 
solar mass BBH inspirals; d) High frequency, which has a 
narrowband tuning at 1 kHz. 

We re-calculate the S/N of the total CBC background 
for H-L using the additional four sensitivity curves, and ob- 
tain 1.18, 0.83, 1,49 and 1.23 for a), b), c) and d) respec- 
tively, in comparison to 1.05 for the standard configuration 
(zero-detuning with high laser power). This shows that mod- 
est improvement of low-frequency sensitivity provide consid- 
erable enhancement in S/N, which is comparable to or even 
greater than that due to the combination of multiple de- 
tector pairs (again for the currently proposed network) . The 
largest value of S/N (for one year observation), which comes 
from the adoption of c), implies that S/N = 3 is achievable 
with an integration time of 4 years in comparison to 8 years 
for the standard option. 

Based on the above analysis, we suggest that the op- 
timal BHBH option offers an appreciable increase in sensi- 
tivities of stochastic searches. To make this suggestion more 
accessible to experimentalists and to extend the above com- 
parison to generic power law GWB models - fiGw(/) — 
fi a (//100Hz) a , we show in Figure the minimum de- 
tectable energy density Q™ n for aLIGO H-L considering 
different tuning options. The values of fi™ n can be easily 
obtained by setting a threshold on S/N and solving the 
equality given by equation (|27fl . We take the integration 
range from 10 Hz to 1 kHz, and consider the range (0 - 



Zero Det, high Power 




1 2 3 4 5 



a 

Figure 9. The minimum detectable energy density 0™ m for 
generic power law GWB models (with indices a) - flGwif) = 
C Q (//100 Hz) a , for one year observation using the aLIGO H-L 
pair. We assume a S/N threshold of 3, and consider five tuning 
options for aLIGO: zero-detuning with high/low laser power, op- 
timized for searches of NSNS/BHBH inspirals and one with high 
frequency narrowband tuning. We refer interested readers to the 
public LIGO document T0900288 for descriptions and technical 
details about aLIGO tunings. 

5) for a. The curves in Figure [9] represent the upper limits 
obtainable by aLIGO, which apply to primordial GWBs (in 
addition to AGWBs), e.g., a = in many early -Universe 
scenarios (see, e.g., Figure 2 in lAbbott et al.ll20091 ). 

5.3 The construction of Qaw(f) templates 

In the previous sections, we have shown that: 

1) for / < 100 Hz, the power law model given by equa- 
tion (|22[1 is a good approximation and requires only three 
parameters. The power law relation holds for three popula- 
tions and thus the total background as well; 

2) above 100 Hz, PN corrections become more notable, and 
different behaviors are expected from other effects such as 
CSFR and mass distributions, making it difficult to predict 
the background spectral properties. 

In this subsection we show that the power law model is suf- 
ficient to be used as search templates for a CBC background 
and is also useful for parameter estimation of the coalescence 
rate and average chirp mass (information other than these 
two quantities can only be extracted from measurements of 
high-frequency peaks). 

Figure [10] illustrates the fractional S/N as a function of 
upper frequency limits for the 10 pairs of advanced detectors 
for the BNS background. We see the S/N has saturated be- 
low 100 Hz due to the suppressing effects of f(f) and the fact 
that the background is "red", i.e., Sh(f) ~ f~ 7 ^ 3 - Quantita- 
tively, a fraction of 99% of the total S/N can be achieved up 
to 51 Hz and (at most) 98 Hz by cross-correlating H-L and 
K-V respectively. Such upper frequencies are slightly higher 
for ET-B (133 Hz) or assuming 7(/) = 1 for aLIGO (128 
Hz) , and could be even lower for ET-D (47 Hz) , as also noted 
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Figure 10. The fractional signal-to-noise ratio (S/N) as a func- 
tion of upper frequency limits for the 10 pairs of advanced de- 
tectors for the BNS background. The "fastest" and "slowest" to 
accumulate 99% of the total S/N are highlighted with thick lines, 
corresponding to cross-correlating aLIGO H-L and KAGRA- 
advanced Virgo (K-V) respectively. 



in iKowalska et aD ((2012) . The exactly same values are ob- 
tained in the cases of BH-NS and BBH due to the similarity 
in fi G w(/) below 200 Hz. 

We then quantify the effectiveness of the proxy of a 
power law model to the CBC background below 100 Hz by 
looking at the following quantity: 



T 



~1 2 {f)S h {f)S' h {f) 
S nl (f)S n2 (f) 



df 



(29) 



which gives the mean value of the cross-correlated signal, 
with A the normalization constant to e nsure (S) = Q a T 
for a power GWB with Qgw(/) = Q / a (|Allen fc Romano! 
1 19991 ). Here S h (f) is the "true" spectral density of the back- 
ground, assumed to be that given by our numerical mod- 
els; S' h (/) corresponds to the template adopted in stochastic 
background searches. A simple power law template results 
in an overestimation of (S) within 2% — 5% for the three 
CBC populations for 10 pairs of advanced detectors. This 
can be further reduced by up to 1% by decreasing the upper 
cutoff frequency from 100 Hz to 50 Hz. 

Overall, we suggest that a simple power law model for 
the CBC background as given by equation (1221) with an ap- 
propriate upper frequency cutoff at 50-100 Hz is sufficient for 
detection and the followed-by parameter estimation of aver- 
age masses and coalescence rates using ground-based inter- 
ferometers. For third-generation detectors like ET, however, 
more accurate models, such as those presented in subsec- 
tions 13.31 and 14.11 will be required to extract information 
such as CSFR, PN effects, and mass distributions. 



6 A FOREGROUND FORMED BY 
SUB-THRESHOLD CBC EVENTS 

When searching for primordial GWBs from the early Uni- 
verse, AGWBs formed by more recent sources could act as 
contaminating foregrounds. One resolution to this problem 
is to subtract individually detected signals from the data. 
This has been demonstrated for the proposed Big Bang Ob- 
server, which has a sufficiently good sensitivity that it can 
resolve and thus subtract away almost all BNS inspirals in 
the U niverse from the overall backg round JCutler fc Harmj 
l2006h . We refer interested readers to lCutler fe Harms! (|2006r ) 
for details of the method and related practical issues, and we 
simple apply this method to ET to estimate the "residual" 
foreground from sub-threshold CBC eventful. 

The optimal method to detect signals with known wave- 
forms is through matched filtering, for which the optimal 
(single-event) sign al-to- noise ratio, p, is given by (see, e.g., 
lAbadie et aflWlOT l: 



p 2 =4 



/ma> 



\Kf)\ 2 

Sn(f) 



df, 



(30) 



where / max is the maximum observed frequency, depending 
on source redshift, component masses and spins (if applica- 
ble). We use the ET antenna pattern function, which goes to 
equation ()26|) for _D e ff and determines the overall amplitudes 
of \h(f)\, for a triangle configuration inc luding three V- 
shap ed detectors (see, e.g., equation (24) of Reg imbau et alj 
2012). Note that the spectral density of the CBC background 
is well below that of the instrumental noise of ET even at 
optimistic rate estimates. Therefore, we do not need to con- 
sider Sh{f) as an additional contribution to S n {f) in equa- 
tion (|30[l . whereas one must do so in the case of Big Bang 
Observer (|Cutler fc Harms! 120061 ). 

We calculate p for each of the simulated CBC events 
in our Monte-Carlo simulation as described in section [4] 
those loudest events resulting in p ^ p t h = 8 are discarded 
("subtracted") to estimate a residual noise. Before moving 
forward to discussions of ET's potential in removing the 
CBC background through a subtraction process, we present 
in Table [4] detection rates (which are conveniently obtained 
in our simulations) of CBC sources for ET. The calcula- 
tions i mprove the ap p roxim ation (for advanced detectors) 
used in lAbadie et al.l (|2010h with the following considera- 
tions: a) cosmic evolution of coalescence rates, the standard 
ACDM cosmology and cosmological redshifts; b) observa- 
tional NS/BH mass distributions; c) complete waveforms 
that include PN amplitude corrections. While these effects 
may not be important for detect ion rate predictions for ad- 
vanced detectors as discussed in lAbadie et al.l (|2010h , they 
must be considered for ET due to its 1000 times larger ac- 
cessible volume. 

Based on results presented in Table [4] we suggest that 
the realistic CBC detection rates for ET be 10 5 (BNS), 10 4 
(BH-NS) and 10 4 (BBH) g iven the current rea listic coales- 
cence rate predictions (see lAbadie et al1l2010l . for details). 
Note that ET will have an overall detection efficiency (de- 
fined as ATdet/iVtot) of ~ 10% (BNS), ~ 40% (BH-NS) and 



As the detection horizon of advanc ed detectors such as aLIGO 
is at most z ~ 0.4 for BBH systems (Abadic et al. 2010), a sub- 
traction process can not reduce the foreground appreciably. 
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Table 4. Detection rates (JVdet) °f CBC sources for ET. 





CSFR 


A/tot 

(yr -1 ) 


G 


P 


(yr" 1 ) 

E 


TG 


BNS 


HB 


255 


30.4 

(31.5) 
18.5 
(19.4) 










RE 


335 










HB 
RE 


7.14 

9.59 


3.20 


3.10 


3.45 


3.42 


BH-NS 


(3.27) 
3.07 


(3.17) 
2.90 


(3.53) 
3.36 


(3.50) 
3.33 




(3.27) 


(2.99) 


(3.45) 


(3.42) 




HB 


1.07 


0.95 


0.93 


0.979 


0.97 


BBH 


(0.96) 


(0.94) 


(0.983) 


(0.98) 


RE 


1.49 


1.22 


1.18 


1.277 


1.26 




(1.23) 


(1.19) 


(1.284) 


(1.27) 



Notes: All values are in 10 4 . Ntot is the to tal event rate up to 
z = 6 or z = 15 for CSFR models o f HB llHopkins fc Beacoml 
l2006f) or RE llRobertson fc EIIidl2012T) respectively. We scale the 
number of events above the detection threshold (p ^ 8) in the 
Monte-Carlo simulation (see section [4} according to A'tot/ATmc to 
obtain the detection rate N^ et . We have considered NS/BH mass 
distributions as described in Table [T]— Gaussian (G), Power law 
(P), Exponential (E) and Two-Gaussian (TG); and adopted ET 
sensitivities of two configurations - ET-B and ET-D (values are 
given in parentheses). ET-D gives slightly higher detection rates 
due to a greater low-frequency (/ < 20 Hz) sensitivity. 

~ 85% (BBH), which is independent of ro and weakly de- 
pendent on coalescence rate evolution and sensitivity models 
as shown in Table [4] 

Figure [TT1 compares the results of flGw(/) calculated 
using equation (|25p without and with the subtraction of in- 
dividually detectable events. We see that ET will be able to 
reduce the CBC background energy densities by a factor of 
about 2, 10 and 200 from the BNS, BH-NS and BBH popu- 
lation respectively through a subtraction scheme. The total 
residual foreground is overwhelmingly due to sub-threshold 
BNS merger events and is insensitive to rate evolutionary 
histories and BH mass distributions. The possibility that ro 
for BBH could be much higher than t he value used here , 
e.g., r = 0.36Mpc- 3 Myr _1 found in iBulik et all l|201ll ). 
does not significantly change the level of such a residual 
foreground. Additionally the contribution from a possible 
population of dynamically formed BBHs to such a residual 
foreground is negligible, as ET will be able to detect these 
sources out to much larger distances than the field popu- 
lation of BBHs consi dered in this work du e to significantly 
higher chirp masses (Sadowski ct al. 20081). 

Figure [12] compares the noise power spectral densities 
of future terrestrial detectors with the spectral densities of 
the CBC residual foreground, the upper bound of th e GWB 
from NS r-mode instabilities (see lZhu et al.ll2011al . for de- 
tails) and a putative primordial GWB fro m the very early 
inflationary Universe (| Abbott et al.l [2~00°l ) . It implies that, 
without considering other types of AGWB sources, while the 
NS r-mode GWB is unlikely to mask the highly-sought in- 
flationary background, the contribution from sub-threshold 
CBC events could be a challenging issue because these sig- 
nals are beyond the capability of current data analysis meth- 




/ (Hz) 

Figure 11. The energy density parameter f2ow(/) of the to- 
tal CBC background and its contributions by the BNS, BH-NS, 
and BBH populations, without and with (labeled as sub) the sub- 
traction of individually detectable sources for ET assuming ET-B 
sensitivity. Note that: a) the total residual "noise" is dominated 
by sub-threshold BNS merger events; b) ET will be able to reduce 
the foreground level by a factor of about 2, 10 and 200 from the 
BNS, BH-NS and BBH population respectively; c)the reduction 
shown here is optimistic because of the assumption that every 
the oretically detectable sin gle signals will be perfectly subtracted 
(see lCutler fc Harms>200r3 . for detailed discussion about potential 
subtraction errors). 

ods and will always add up to act as an additional "noise" 
component in the data. 



7 TIME-FREQUENCY PROPERTIES 

The analysis in the previous section is based on the expec- 
tation that individual CBC events contributing to a back- 
ground have different amplitudes in the data due to a dis- 
tribution over source distances, orientations and sky posi- 
tions. Rigorously one need to track the number of sources 
contributing in the relevant frequency intervals. For CBC 
events, the duty cycle is frequency dependent - individual 
signals stay much longer at low frequencies, leading to much 
smaller duty cycle for incr easing frequencies, as pointed out 
in recent studies (see, e.g.. lRosadoll201ll ; IWu et al.ll2012l ). 

Practically we want to know the critical frequency, / c , 
above which individual signals do not simultaneously occupy 
the same frequency interval. Without going into specific de- 
tails of the calculations, we provide the following relation: 

/ f c V 1/3 _(Af\ f r \ / Mc V 5/3 

1,15 Hz,/ V 1Hz / \lMpc- 3 Myr-V \ 1 Mq J 

(31) 

where A/ is the size of frequency interval relevant to the 
analysis. We consider a reference value of 1 Hz for A/, while 
in practice it could be as small as the frequency resolution of 
an experiment. For ground-based interferometers such a res- 
olution is given by 1/ AT, where AT is the time duration of 
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Figure 12. The spectral densities, Sh(f), of the residual fore- 
ground formed by sub-threshold CBC events, the GW B from NS 
r-mode instabilities (shown here is the upper bound, see lZhu et alj 
l|2011ah for details) and a putative primordial GWB from infla- 
tion (with a constant energy density parameter Hq^v = 10 -14 ). 
All are compared against noise power spectral densities, S n (f), 
of second (aLIGO, K - KAGRA, V - advanced Virgo) and third 
generation (two possible configurations of the planned Einstein 
Telescope, ET-B and ET-D) ground-based GW interferometers. 
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Figure 13. The spectrograms for a simulated BNS background 
signal (upper panel) and for the same background signal plus 
an arbitrary amount of Gaussian white noise (lower panel). In 
the simulation we increase the coalescence rate by a factor of 
50 and scale down the final amplitude by the same factor. The 
colors or intensities in the spectrogram effectively tells the relative 
signal/noise power spectral densities at given time and frequency 
instances. The lower panel is only for illustration and should not 
be related to actual detection prospects. 



short data segments which are used in cros s correlation anal- 
ysis and is typically of o rders of seconds l|Allen fc Romano] 
Il999l ; lAbadie et afll2012l ). In equation fl3TJ we neglect the 
effect of mass distribution without loss of generality. 

Once we know the critical frequency f c , the duty cycle 
function £(/, A/) is exclusively determined through: 



e(/,A/) 



(Al 
V 1Hz 



-11/3 



(32) 



Note that: a) the duty cycle fun ction is comparable to the 
overla p function in iRosadd ([20111 ) and the duty cycle param- 
eter in lWu et all (|2012l ): b) the contribution of post-inspiral 
emission to £(/, A/) is negligible due to much shorter dura- 
tions; c) the above two convenient relations are applicable 
to three CBC populations up to a few hundreds Hz. 

The physical meaning of the duty cycle function is the 
(statistically average) number of intersections of the tracks 
of individual inspirals in time-frequency plane at a given fre- 
quency interval. We demonstrate this in the upper panel of 
Figure [13] by plotting the spectrogram of a simulated time 
series of BNS inspirals up to z — 6. This allows one to 
visualize the unique time-frequency properties of the BNS 
background: a) above a few tens Hz, individual chirps are 
only occasionally present, separated in both time and fre- 
quency domain; b) a large number of overlapping signals be- 
low ~ 20 Hz create a continuous and stochastic background, 
with the colors showing the "redness" of the background. 
The simulation f ollows the same procedure as described in 
iRegimbau et al l (|2012h except that we increase ro by a fac- 
tor of 50 and scale down the final amplitude by the same 
factor. 

For illustration we add an arbitrary amount of Gaussian 



white noise to the simulated background, and the spectro- 
gram is shown in the lower panel of Figure [13] As most of 
the signals are deeply buried in detector noise, one can only 
resolve the strongest events which are above the detection 
threshold, just like a few "chirping" structures apparent in 
the noisy spectrogram. Note that the comparison of signals 
and noise in this plot are not representative of the actual 
detection prospects since the added noise is not comparable 
to (and obviously weaker than) the instrumental noise of 
current and future ground-based detectors. 

The conclusion from equation (|31|l and Figure[T3]is that 
at each 1 Hz frequency bin above 15 Hz, we will be observ- 
ing individual BNS inspirals rather than a background. The 
equivalent critical frequency is of order 4 (2) Hz for the BH- 
NS (BBH) population. Furthermore, when an array of de- 
tectors with moderate angular resolution is used to observe 
the CBC background, the background must be considered 
as a set of discrete transient sources randomly distributed 
across the sky. However, this does not necessarily mean the 
GW emission due to various CBC populations as a whole 
can not be detected as a background by the standard cross 
correlation method (as we have pointed out in the beginning 
of subsection 15. 1|1 . In practice, the underlying numerous in- 
dividual signals are averaged when cross-correlating data of 
length from months to years, and it is possible to recover the 
theoretical expected spectral density as de monstrated in a 
recen t mock data challenge study for ET (fRcgimba u et al] 
120121) . In the same way, sub-threshold transient signals re- 
main in the data as an additional noise component which 
could obscure the primordial GWBs and other AGWBs. 



GWB from compact binary coalescence 17 



8 CONCLUSIONS 

In this paper we first reviewed the formalism of the calcu- 
lation of £2gw(/) and developed a practical model for AG- 
WBs - a power law energy spectrum dEGw/df ~ nat- 
urally leads to Hgw(/) = ^a.f a where £l a depends almost 
exclusively on the local rate density ro and total amount 
of radiated GW energy A_Egw- Such a model allows one 
to quickly evaluate uncertainties in estimates of the back- 
ground strength and the associated detectability. 

We have provided updated estimates of the spectral 
properties of the CBC background formed by populations of 
BNS, BH-NS and BBH systems. By systematically investi- 
gating effects of CSFRs, delay times, NS/BH mass distribu- 
tions, and using up-to-date analytical complete waveforms 
including PN amplitude corrections, we showed that: 

1) Effects of CSFRs and delay times are linear below 100 Hz 
and can be represented by a single parameter J2/3 with an 
uncertainty ~ 2; 

2) PN effects cause a small reduction of 0.gw (/) from a f 2 ' 3 
power law function above a few tens Hz; 

3) Below 100 Hz, SlGw(f) can be approximated by a f 2 ^ 3 
power law function, with the magnitude determined by only 
three parameters - the local coalescence rates ro, the aver- 
age system mass (M^ 3 ) plus J2/3. In particular, within this 
frequency range f2ow(/) does not depend on chirp mass dis- 
tributions. This finding, which was also obtained indepen- 
dently i n a recent study usin g a population-synthesis ap- 
proach (|Kowalska et alj 12012). is important for parameter 
estimation of this background; 

4) A variety of features at high frequencies (> 200 Hz), 
e.g., different peak frequencies and widths of S1gw(/), are 
expected from different CSFRs, delay times, and mass dis- 
tributions. Measurements of the peaks will be rewarding al- 
though challenging due to the small contribution (less than 
1%) to S/N by the high frequency signal; 

5) The post-merger emission of BNS coalescences could 
considerably enhance the peak of the BNS background at 
around 1-2 kHz, but will not alter the background spectrum 
below 300 Hz. While this contribution to a GWB may be 
too weak to be detectable even for ET, the latter fact is 
advantageous for parameter (ro and (M^ 3 )) estimation by 
measuring only the low-frequency power law spectrum. 

Using updated estimates of Q,Gw(f), we revisited the 
issue on the detectability of this background signal. Assum- 
ing a detection target of the total background contributed 
by three CBC populations for a worldwide network of ad- 
vanced detectors, we showed in Figure [S]the accessible "rate 
space" of the local coalescence rates ro (in Mpc -3 Myr" 1 , 
with the value of BH-NS fixed at 0.03), implying: 

1) A combination of a BNS population at the realistic rate 
of ro = 1 and a BBH population at a rate of ro = 0.1 will 
give rise to a detectable background signal; 

2) Either a BNS rate of r = 2.7 or a BBH r = 0.16 will 
be necessary for detection, when BBH or BNS has very low 
coalescence rate (note that the chosen ro for BH-NS ensures 
a negligible contribution). 

In both cases, recent optimistic rate estimates for BBHs pro- 
vide interesting detection prospects for a CBC background. 
The above quoted values are for optimally combining a net- 
work of 5 advanced detectors. Such an optimal combination 
gives 30% improvement in detectability over aLIGO H-L. 



This is way below the common expectation l|Wu et alj|2012l : 
iKowalska et all l2012l l that such a network could perform 
as well as two co-located and co-aligned aLIGO detectors, 
which gives a 3-fold improvement on H-L. In the l atter case 
our re sults are consistent with those presented in IWu" et all 
|2012l ). showing it is likely that at the realistic rate a BNS 
background may be detected within one year observation 
using two co-located aLIGO interferometers. 

We emphasize that the somewhat "disappointing" per- 
formance of a network of detectors is due to effects of the 
overlap reduction functions for the current configurations 
of the advanced detector array - the large separations be- 
tween pairs of detectors, of the orders of 10 4 km (except H-L, 
3000 km), result in very modest correlation of background 
signals above 20 Hz (50 Hz for H-L). This further implies 
that stochastic background searches can benefit significantly 
from a pair of closely spaced detectors, with separation cho- 
sen to be both within one reduced wavelength (about 300 
km for 150 Hz) and relatively large to ensure that their noise 
sources are largely uncorrelated. 

We found that 99% of the S/N can be obtained by con- 
sidering only the contribution up to 50 Hz (aLIGO H-L) or 
at most 100 Hz (KAGRA-advanced Virgo). Two main im- 
plications for advanced detectors are: 

1) Only the low frequency part is important for detection; 

2) Improvement on the sensitivity below 50 Hz is beneficial 
for detection. 

We conclude that a simple power law model as given by 
equation (|22|) with an upper frequency cutoff of 50-100 Hz 
is sufficient for background searches. Since the model is gen- 
eralized to three CBC populations and only requires three 
parameters, it could prove useful to constrain or estimate 
these parameters with future stochastic searches - partic- 
ularly one can marginalize over a uniform distributed J 2 /z 
to obtain confidence levels of ro and {M^ 3 }. In addition, 
our generalized model can also be used to identify the rela- 
tive contribution from different populations from the likely 
detection of a CBC background. This will further require 
combination of stochastic ba ckground measurem ents with 
CBC single event detections (Ma ndic et alj|2~012l ). Regard- 
ing the above point 2), we specifically showed that for the 
CBC background the aLIGO tuning configuration offering 
the best low-frequency sensitivity (which is optimized for 
BBH inspiral searches) will provide a 50% enhancement in 
the achievable S/N against the standard sensitivity (zero 
detuning with high laser power), and such an improvement 
is even better than that due to an optimal combination of 
the currently proposed detector network (which comprises 5 
advanced detectors). We further compared the sensitivities 
of stochastic searches using different aLIGO tuning options 
to generic power law GWB models in Figure [JJ] The results 
show that aLIGO H-L will be able to detect a GWB with a 
S/N above 3 for fio ^ 1-87 x 10 -9 (assuming a flat energy 
spectrum) with one year observation at the standard sensi- 
tivity, and this limit could be reduced down to 1.24 x 10 -9 
using the optimal BHBH option. 

For third generation detectors like ET, the background 
will be easily detectable, with a S/N from tens up to 
hundreds contributed by individual populations. The high 
achievable S/N will open up the possibilities to: a) enable 
different populations to be disentangled; b) measure PN ef- 
fects; c) probe mass distributions and rate evolutionary his- 
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tories by measuring the (high-frequency) peaks of the back- 
ground energy spectra. To gain more insights about how 
these information can be extracted from background mea- 
surements, models presented in this study can be further 
improved in the following ways: 

1) Contribution from possible populations of dynamically 

formed BBHs and/or binaries involved with one or two 

l '. 1 1 1 

intermediate-mass BHs (see Ab adie et aL||2010| , and refer- 
ences therein for details) should be considered. Due to sig- 
nificantly higher masses of such systems, their contribution 
could peak at a few tens Hz and potentially distort the power 
law relation for the three normal CBC populations; 

2) More accurate complete waveforms are required. In this 
regard, the three types of BH- NS waveforms corresponding 
to different merger processes l|Shibata fc Taniguchil [20111 ) 
are of particular interest and can be used to investigate how 
the information of NS equation of state is encoded in the 
background signal. 

We demonstrated that ET will be able to reduce the 
contributions to a GWB from the BNS, BH-NS and BBH 
populations respectively by a factor of 2, 10 and 200 through 
the subtraction of individually detectable events but there 
is a strong residual foreground dominated by sub-threshold 
BNS merger events. Such a foreground, at the level of 
flaw ~ 10 -10 in the (1-500) Hz frequency range, can hardly 
be removed and must be considered in future terrestrial 
searches of primordial GWBs and other AGWBs. 

We finally discussed the unique properties of the CBC 
background - well defined continuously rising tones, local- 
ized directions and well defined average spectral density. 
These have not so far been fully exploited by stochastic 
background searches. We believe that new algorithms could 
exploit these properties to go beyond the standard cross 
correlation limit that applies only to true stochastic back- 
grounds. 
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